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Abstract 

We consider an ultracold quantum degenerate gas in an optical lattice inside a cavity. This 
system represents a simple but key model for "quantum optics with quantum gases," where a 
quantum description of both light and atomic motion is equally important. Due to the dynamical 
entanglement of atomic motion and light, the measurement of light affects the many-body atomic 
state as well. The conditional atomic dynamics can be described using the Quantum Monte Carlo 
Wave Function Simulation method. In this paper, we emphasize how this usually complicated 
numerical procedure can be reduced to an analytical solution after some assumptions and approx- 
imations valid for macroscopic Bose-Einstein condensates (BEC) with large atom numbers. The 
theory can be applied for lattices with both low filling factors (e.g. one atom per lattice site in 
average) and very high filling factors (e.g. a BEC in a double- well potential). The purity of the 
resulting multipartite entangled atomic state is analyzed. 
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I. INTRODUCTION 



Since the first observation of the Bose-Einstein condensation (BEC) in 1995, physics of 
ultracold quantum gases has become a well established field considering various collective 
quantum states of bosonic and fermionic atoms trapped in optical potentials. However, in 
the majority of both theoretical and experimental works, the role of light is reduced to a 
classical axillary tool for creating fascinating quantum atomic states. On the other hand, 
recent experimental achievements, where a quantum gas was loaded into a typical quantum 
optical setup (high-Q cavity) [l, 2, ^,0,0], provides a challenge to develop a theory of novel 
phenomena, where the quantizations of both light and atomic motion pla y equall y important 
roles. Recently, we have contributed to this filed 
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At present, the main procedure to measure the properties of ultracold atoms is the time- 
of-flight method, where the trapping potential is switched off and the interference of matter 
waves is observed. This method is completely destructive as the preparation of a new atomic 
system is necessary for each measurement. As we have shown in 7J, |8|, |9J, the light scattering 
provides a much less destructive method to measure atomic properties, i. e., the quantum 
nondemolition (QND) measurement scheme. Light scattering does not destroy the atomic 
system, hence, many consecutive measurements of light can be done using the same atomic 
sample without preparing a new one. 

However, as Quantum Mechanics states, any measurement affects the quantum state. 
So, even if the QND measurement does not destroy the atomic system, it still changes the 
quantum states of atoms. In Refs. 

BBS 

, ll3j we presented some relations between the 



expectation values of the atomic and light variables. Essentially, one needs many measure- 
ments to obtain the expectation value of some quantum variable. Thus, even for optical 
measurements, one needs to prepare the initial state several times (in contrast to the time- 
of-flight schemes, the initial state can be prepared with the same atoms). 

In Refs. UlLdJ, we put a different question: the continuous measurement of light scat- 
tered from the atoms was considered. Due to the light-atom entanglement, the measurement 
of light affects the many-body atomic state as well. Thus, the quantum back-action of light 
measurement on the atomic state was analyzed. The conditional dynamics of the atomic 
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FIG. 1: Setup. A lattice is illuminated by the transverse probe ao and probe through a mirror 
■q. The photodetector measures photons of the cavity mode a\ leaking out of the cavity. Due to 
the quantum back-action, the light measurement leads to the modification of the atomic quantum 
state. 

state due to the light measurement was presented. In contrast to calculating the expectation 
values, such an approach gives the evolution of a system at a single quantum trajectory. This 
evolution should be first seen by experimentalists before they do multiple measurements and 
average them. Moreover, this measurement can be used as a method to prepare particu- 
ar atomic multipartite entangled states thanks to the quantum measurement back-action 

A standard procedure to analyze the conditional evolution of the atomic system observing 
light is the Quantum Monte Carlo Wave Function (QMCWF) simulation method. Usually, 
it is a basis for numerical simulations of the quantum dynamics in open dissipative systems. 
In this paper, we show how this usually complicated numerical procedure can be reduced to a 
simple analytical solution after some assumptions and approximations valid for macroscopic 
Bose-Einstein condensates (BEC) with large atom numbers. The purity of the resulting 
multipartite entangled atomic state is analyzed. 
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II. THEORETICAL MODEL AND QUANTUM MEASUREMENTS 



We consider the model presented in Refs. [Ill, ll2| . Using the approximation of macro- 
scopic atomic ensemble with the large atom number, which is relevant for present experi- 
ments, will enable us to obtain a simple analytical solution in the next section. 

We consider (cf. Fig. 1) N ultracold atoms in an optical lattice of M sites formed by 
strong off-resonant laser beams. A region of K < M sites is also illuminated by a weak 
external probe, which is scattered into a cavity. Alternatively, this region is illuminated by 
the cavity field appearing due to the presence of the probe through the cavity mirror. 

We use the open system approach for counting photons leaking the cavity of decay rate 
k. When a photon is detected, the jump operator (the cavity photon annihilation operator 
ai) is applied to the quantum state: \^ c (t)) — > ai\^ c (t)). Between the counts, the system 
evolves with a non-Hermitian Hamiltonian. Such an evolution gives a quantum trajectory 
for |\l/ c (£)) conditioned on the detection of photons. 

The expression for the initial motional state of atoms reads 

i*(°)> = ••,?*>. C 1 ) 

q 

which is a superposition of Fock states reflecting all possible classical configurations q = 
{?!)••> <?m} °f N atoms at M sites, where qj is the atom number at the site j. As we 



have shown in Refs. 
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12| . the solution for conditional wave function takes physically 
transparent form, if the following approximations are used: atomic tunneling is much slower 
than light dynamics, the probe waves are in the coherent state, the first photon is detected 
after the time 1/k. The conditional state after the time t and m photocounts is given by 
the quantum superposition of solutions corresponding to the atomic Fock states in Eq. (JT]): 

l*e(m,f)> = ^yE< e * ,Wc Sl^'-'^)l^)' ( 2 ) 



7] - iU l0 a D q 
aq iiUnD^-AJ+K' [) 
$ 5 (t) = -\a q \ 2 nt + (r]a* q - iU 10 a D q w a* q - c.c.)t/2, (4) 

where a q is the cavity light amplitude corresponding to the classical configuration q. It is 
simply given by the Lorentz function ([3]) well-known from classical optics, where a® is the 
external probe amplitude, rj is the amplitude of the probe through a mirror; Ui m = gig m / A a 



(l,m = 0, 1), where g^o are the atom-light coupling constants, A a = oj\ — u a is the cavity- 
atom detuning; A p = u p — uj\ is the probe-cavity detuning. D\ m = Ylf=i u *{ r j) u m{Yj)<lj are 
the probe-cavity coupling coefficient and dispersive frequency shift that sums contributions 
from all illuminated atoms with prefactors given by the light mode functions Uo,i(r). Except 
the prefactors associated with m photodetections a™, the components of quantum super- 
position in Eq. (j5J) acquires the phases contained in $ g (t) (jlj). F(t) is the normalization 
coefficient. 

For several particular cases, the solution for the time-dependent probability distribution 
of atoms corresponding to the state (J2J) can be simplified further [ll|, Q. If the probe, 
cavity, and lattice satisfy the condition of the diffraction maximum for light scattering, the 
probability to find the atom number < z < N in the lattice region of K sites is given by 

p(z im ,t)=z 2m e-^p (z)/F 2 , (5) 
F 2 = J2 z2m e~ z2r Po(z), 

z 

with r = 2\C\ 2 Kt, C = iUi a Q / (iA p — k), po(z) is the initial distribution, and F provides the 
normalization. The light amplitude corresponding to the atom number z is a z = Cz. 

If the condition of a diffraction minimum is satisfied, the probability to find the atom 
number difference between the odd and even sites in the lattice region of K = M sites is given 
by the same Eq. (jSJ), but with a different meaning of the statistical variable —N < z < N. 

When the time progresses, both m and t increase with an essentially probabilistic relation 
between them. The Quantum Monte Carlo method, which establishes such a relation thus 
giving a quantum trajectory, consists in the following. The evolution is split into small 
time intervals 5tt. In each time step, the conditional photon number is calculated using the 
probability distribution flS}: 

(a\ ai ) c (t) = \C\ 2 z Mz, ™, t), (6) 

z 

which is proportional to the second moment of p(z,m,t). The probability of the next, 
(m + l)th, photocount within this time interval P m +i — 2n{a\ai) c Sti is then compared with 
a random number < q < 1 generated in advance, thus, deciding whether the detection (if 
2n{a\ai) c 8ti > e<) or no-count process (otherwise) has happened. 
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III. MACROSCOPIC QUANTUM GASES 



In this section, we consider a case, where the initial atomic state is a macroscopic super- 
fluid (SF) with the atom number N ^> 1. Note, that the total number of lattice sites M 
and the number of sites illuminated K can be any. Thus, the theory presented below can be 
applied for lattices with both low filling factors (e.g. one atom per lattice site in average) 
and very high filling factors (e.g. a BEC in a double- well potential). 
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We will start with the case of a diffraction minimum. As shown in Refs. 
probability function (jSJ) to find the atom number difference z between the odd and even sites 
shrinks to a doublet with the peaks at z\$ = y^m/r, which corresponds to the generation 
of the Schrodinger cat state. 

For the SF state the probability to find the atom number at odd (or even) sites z [z = 
(z + N)/2 as the atom number difference is z and the total atom number is N] is given by 
the binomial distribution 

PsfCz) = wrhy. (1) i 1 -^) ' (7) 

where Q is the number of odd (or even) sites. For even M, Q = M/2 and Eq. (J7]) simpli- 
fies. For a lattice with the large atom number N ^> 1, this binomial distribution can be 
approximated by a Gaussian function. Changing the variable as z = 2z — N we obtain the 
Gaussian function for the probability to find the atom number difference z: 

Psf{z) = ~j=^z ^ (8) 
V zira 



with the zero mean z and a = yN giving the full width at half maximum (FWHM) 



2oV2 In 2. The variance of the atom number difference in the SF state is o = N. 

For the large atom number, the summations in Eqs. (jSJ) and (J6j) can be replaced by the 
integrals over all — oo < z < oo, which gives the following probability of the next photocount: 
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Pm+i = — ^ _ z2 $n, (9) 

f z 2m e~ z2r e ^dz 

J — oo 



where 6ri = 2\C\ 2 K,5ti. Taking into account the following relation 32J 
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the integrals can be calculated and the probability of the next photocount reads 

Thus, we see that the Quantum Monte Carlo method, which is usually expected to be 
a basis for hard numerical simulations, has reduced to an extremely simple form. After 
splitting our time axis into intervals St^ and generating the random numbers < < 1, 
one has simply to substitute the current time r and the photocount number m in the trivial 
algebraic expression (TTTj) and realize if the next photocount happened or not. Proceeding 
this way one establishes the relation between the photocount number and time m(r) at the 
quantum trajectory corresponding to the generated set of random numbers 6j. Knowing the 
relation between m and r, one can calculate various conditional expectation values using 
Eq. fflQj) and the complementary expression 



n 



o 



x 2n+1 e~ px dx = ——. (12) 
2p n+i v ) 



We now switch to the case of the light detection at the direction of a diffraction maximum. 



As shown in Refs. 
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121 ] , in this case one measures the atom number z in the lattice region 
of K illuminated sites. 

For the initial SF state the probability to find the atom number z at the lattice region of 
K sites is given by the binomial distribution 

For a lattice with the large atom number N 3> 1, it can be approximated as a Gaussian 
distribution 

1 (z-zq) 2 

Psf(z) = -7=^e (14) 
VZira 



with the mean atom number z = NK/M and a = y/N(K/M)(l - K/M) giving FWHM 
2<7 V2 In 2. The atom number variance in the SF state is a 2 . 

Similarly to the case of a diffraction minimum, for the large atom number, the summations 
in Eqs. (jSJ) and ([6]) can be replaced by the integrals over — oo < z < oo, which gives the 
following probability of the next photocount: 

L 00 z 2m+2 e-^e~^dz , 

P m+ i = ^ — g Sn. (15) 

J °° z 2m e~ z2r e 2 CT ^ dz 



Taking into account the following relation 321 ]: 



^e-^ dx = ^r_^_ [q ^ ] (16) 



1 f¥ d n ~ ] 
~ l p Y p dq 

n E(n/2) 1 / \k 



21 TT f q\ v-^ 1 / p 

me p 1 



p \P J — 2/c)!/c! \4g 2 

where E(n/2) is the integer part of n/2, the integrals can be calculated and the probability 
of the next photocount reads 

m+l m , 

P m+1 = (2,n + l)(2 m + V £ (2m + 2 _ afc)|t , / £ (2m _ afc)!t! . < 17 > 
where the parameters are 

z 1 



2a 2 (r + l/(2(7 2 )) 2(1 - K/M)(r + 1 /(2a 2 )) ' 
u r+l/(2a 2 ) 4 iO 2 



Expression (|17[) is more complicated than Eq. (fTTj) for the diffraction minimum. However, 
it is also very simple as it includes only summation over the photocount number. Thus, 
we were able to replace the summation over the atom number, which can rich the values of 
N = 10 5 , and even the numerical integration over the atom number. For the far off-resonant 
interaction considered here, the number of photocounts m will be many orders of magnitude 
less than the atom number. So, the sum in Eq. ( fTTl) will contain only a small number of 
terms. 



IV. PURITY OF THE FINAL STATE 



As shown in Ref. 121 ] . under some conditions, the entangled light-matter state can 



collapse to a macroscopic superposition state, which has a form 

|*c> = ^[e fr k>K> + e-*|*2>K>], (18) 

where 7 is some phase, | ^1,2) are the atomic Fock states with precisely known atom num- 
bers Z\ and Z2- \ot Zl 2 ) are the corresponding coherent states of light such that the light 
amplitudes have equal absolute values, but opposite phases: a Zl = \a Zl | exp {itp) and 
a Z2 = \a Zl \exp(-i(p). 



Expression (fTB|) gives a macroscopic superposition state, where light and atoms are en- 
tangled. One possibility to disentangle them is to switch off the probe and count the photons 
leaking out of the cavity. In this case both \a Zl ) and \a Z2 ) will approach the vacuum state 
|0) and the light and matter will factorize. The difficulty of such a method consists in the 
necessity to count almost all photons to get finally an atomic state of high purity [121 ] . In 
this paper, we analyze the state obtained by tracing out the light field directly in Eq. ({TBI) , 
without waiting for photon leakage. This corresponds to a scheme, where the measurement 
of light is not performed. 

The trace over the light variables can be calculated using the photon Fock basis, in which 
the coherent light states read 



12/, A \a z A n e± n * 



1^21,2/ 



) = e -l^l 2 /2^ l^l 1_ \n>. (19) 



v n 

n=0 



The density matrix of the state after tracing out the light is p = J2^=o( n \^c) (^c\ n )- Using 
the expression 



\2n 



P«!_ e -|<^IV*HP = e l«. 1 l 2 (e- 2 --l) ) (20 ) 
n=0 ' 

the density matrix takes the form 

p= l - (\ Zl ){z x \ + \z 2 )(z 2 \ + e\ a ^ e ^-^e 2 ^\z 1 )(z 2 \ + e^il 2 ^- 1 ^- 2 ^)^) (21) 

The quantity characterizing how close is a mixture state to a pure state is the so-called 
purity: P = Tr(p 2 ). For a pure state it is maximal and equal to 1, while for a maximally 
mixed state it is minimal and equal to 1/2 (in our case of the two-component states). The 
purity of the state (1211) is given by 

P= ~ (i + e-^lW^ _ (22) 

The purity depends on the amplitude of the coherent light states in Eq. ffT8l) and the phase 
difference between them. In a trivial case, where two coherent states are indistinguishable 
(ip = 0), the purity is maximal, p = 1, and the state is pure (however, in this case, z\ = z 2 and 
the state is not a macroscopically entangled one). In non-trivial cases, where the coherent 
states differ by the phase 2<p, the purity decreases with increase of the light amplitude and 
phase difference between them. One can estimate the minimal possible purity as follows. In 
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a coherent state, the uncertainty of the X-quadrature is 1/2. Thus, two coherent states can 
be well distinguished, if \a Zl \ sinyj > 1/4. Substituting the minimal value 1/4 in Eq. fl22l) . 
one sees that the maximal purity can reach 0.89, which is a rather high value. 



V. CONCLUSIONS 



We have considered the quantum measurement of light scattered from an ultracold quan- 
tum gas. Due to the light-matter entanglement, the back-action of light measurement mod- 
ifies the atomic state as well. We applied our theory developed in Refs. [ll|, Q for the 
case of a macroscopic BEC with a large atom number. We have shown that the Quantum 
Monte Carlo simulation method gets very simple in that case and the calculations can be 
carried out even analytically. The purity of the macroscopic superposition atomic state was 
analyzed and shown to be able to reach the high values. 
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